Imperial/TP/99-00/?? 
hep-ph/9912?? 

ETEX-ed on February 1, 2008 



^ Fluctuations and Phase Transition 

o; Dynamics l\ 

(N 

^ R. J. Rivers 

> . 

, Blackett Laboratory, Imperial College, London SW7 2BZ 

0\ 

o 
o 
o 

> 
X 



Abstract 

Kibble and Zurek have provided a unifying causal picture for the appearance of 
classical defects like cosmic strings or vortices at the onset of phase transitions in 
relativistic QFT and condensed matter systems respectively. In condensed matter 
the predictions are partially supported by agreement with experiments in superfluid 
helium. We provide an alternative picture for the initial appearance of defects that 
supports the experimental evidence. When the original predictions fail, this is 
understood, in part, as a consequence of thermal fluctuations (noise), which play 
a comparable role in both condensed matter and QFT. 



^To be published in the proceedings of the 4th Conference on Quantum and stochastic gravity. String 
cosmology and Inflation, Peyresq, September 1999 
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1 Overview 



In this talk I want to consider the emergence of 'classical' field configurations - topological 
defects - after a phase transition, and the extent to which thermal fiuctuations can inhibit 
this process. This is of particular interest in the early universe, for which we expect a 
sequence of transitions from a very symmetric initial state, and in which the presence of 
classical defects can have important astrophysical consequences. 

The relevance of topological defects is that when symmetry breaks, it does not do so 
uniformly. At the very least, the field is uncorrelated on the scale of the causal horizon 
at any time. Since a broken symmetry is, necessarily, characterised by degenerate vacua, 
the choice of different vacua in domains in which the fields are uncorrelated will lead 
naturally to topological defects between them as the field does its best to order on large 
scales. The nature of the defects depends on the relevant homotopy group of the ground- 
state manifold. The most acceptable defect on cosmological grounds is the 'cosmic string' 
- a generalised field vortex - which may have played in role in structure formation. 

Given a theory that permits vortices, at some time after the transition we expect to 
find a network of them, behaving essentially classically as Nambu-Goto strings, intersect- 
ing, chopping off loops which decay, and straightening segments to reduce field gradients. 
Simple calculations suggest late-time scaling solutions, with similarity on a wide range 
of scales. 

The details do not concern us here. What interests us is how this collection of essen- 
tially classical objects, which can be observed directly, in principle, came into existence. 
The simplest question, that we shall address here, is what is the density of cosmic strings 
(or other defects) at the time of their appearance? 

The early universe is very hot, but such a problem requires us to go beyond equilibrium 
Thermal Field Theory. In practice, we often know remarkably little about the dynamics 
of thermal systems. For simplicity, I shall assume scalar field order parameters, with 
continuous transitions. In principle, the field correlation length diverges at a continuous 
transition. In practice, it does not since there is not enough time. One possibility is that 
the separation of 'defects' is characterised by the correlation length when it checks its 
growth. If this were simply so, a measurement of defect densities would be a measurement 
of correlation lengths. Estimates of this early field ordering and its contingent defects in 
the early universe have been made by Kibble[|l], ||, using simple thermal[jl| arguments or 
causal arguments different from the one above (although that is also due to Kibble0). 

There are great difficulties in converting such predictions for the early universe into 
experimental observations since, but for a possible stray monopole, we have no direct 
evidence for them having existed |. Zurek suggested 0] that similar arguments to those 
in were applicable to condensed matter systems for which direct experiments on defect 
densities could be performed. This has lead to considerable activity from theorists work- 
ing on the boundary between QFT and Condensed Matter theory and from condensed 

^Although this does not impede our ability to make predictions for defect-driven fluctuations in the 
CMB, for example 
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matter experimentalists. To date almost all experiments have involved superfiuids, for 
which vortices can be produced readily. All but one experiment is in agreement with 
these simple causal predictions and we shall pay particular attention to this one failure 
of prediction. In this talk I shall 

• review the Kibble/Zurek causality predictions for initial correlation lengths and 
defect densities. 

• summarise the results of the condensed matter experiments and present an alter- 
native picture for the onset of defect production for condensed matter systems. I 
shall then show how this alternative picture gives essentially the same results as the 
Zurek picture for those condensed matter systems for which there is experimental 
agreement. 

• provide an explanation for why some condensed matter experiments will be in 
disagreement with Zurek's predictions, including the experiment in question. We 
shall suggest that the prediction fails, in part, because of the presence of thermal 
noise. 

• use these ideas to address the more comphcated problem of the appearance of 
'classical' defect configurations in QFT in the light of Kibble's predictions, and the 
role of thermal noise in them. 

2 When symmetry breaks, how big are the smallest 
identifiable pieces? 

Defects in the large-scale ordering of the field can only appear once the transition has 
taken place. If it is the case that defect density can be identified simply from the field cor- 
relation length, the maximum density (an experimental observable in condensed matter 
systems, although not for the early universe) will be associated with the smallest identi- 
fiable correlation length in the broken phase once the transition has been effected. This 
provides the initial condition for the evolution of field ordering. From this viewpoint, we 
can observe the defects by default merely by determining the correlation length at that 
time. 

In order to see how to identify these 'smallest pieces'^ it is sufficient to consider the 
simplest theory with vortices, that of a single relativistic complex scalar field in three 
spatial dimensions, undergoing a temperature quench. In the first instance we assume 
that the qualitative dynamics of the transition are conditioned by the field's equilibrium 
free energy, of the form 

F(T) = I d'x (|V<^r + m'{TM' + Al^r) (2.1) 
■^The title of this section is essentially that posed in recent papers by Zurekp. 
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Prior to the transition, at temperature T > Tc, the critical temperature, m{T) > 
plays the role of an effective 'plasma' mass due to the interactions of (p with the heat 
bath, which includes its own particles. After the transition, when T is effectively zero, 
m^(0) = — < enforces the U{1) symmetry-breaking, with field expectation values 
= 7], 7f = M^/A. The change in temperature that leads to the change in the sign 
of m? is most simply understood as a consequence of the system expanding. Thus, in 
the early universe, a weakly interacting relativistic plasma at temperature T 3> M has 
an entropy density s oc T^. As long as thermal equilibrium can be maintained, constant 
entropy S per comoving volume, S oc sa(t)^, gives T oc a{t)~^ and falling, for increasing 
scale factor a{t). Models that attempt to take inflation into account, however, lead to 
'preheating' that is not BoltzmannianQ. Nonetheless, even in such cases it is possible to 
isolate an effective temperature for long-wavelength modes. This is all that is necessary, 
but is too sophisticated for the simple scenarios that we shall present here. We shall not 
even include a metric in Eq. |2.1| . 

The minima of the final potential of Eq. |2.1| now constitute the circle = r/e*". When 
the transition starts begins to fall into the valley of the potential, choosing a random 
phase. This randomly chosen phase can vary from point to point subject to continuity. 
At late times the failure of the field to be uniform in phase on large scales will lead 
to it twisting around classical 'defects' - solutions to 6F/6(f) = that locally minimise 
the energy stored in field gradients and potentials. Those of interest to us are vortices, 
tubes of 'false' vacuum 0^0, around which the field phase changes by ±27r. In an early 
universe context these are the simplest possible 'cosmic strings'. 

How this collapse takes place determines the size of the first identifiable domains. 
It was suggested by Kibble and Zurek that this size is essentially the equilibrium field 
correlation length C,eq at some appropriate temperature close to the transition. I shall 
argue later that this is too simple but, nonetheless, it is a plausible starting point. Two 
very different mechanisms have been proposed for estimating this size. 

2.1 Thermal activation 

In the early work on the cosmic string scenario an alternative possibility to simple causal- 
ity was to assume [|I| that initial domain size was fixed in the Ginzburg regime by the 
correlation length at that time, rather than the causal radius. By this we mean the 
following. Once we are below Tc, and the central hump in V{(f)) = m^(T)|(/)p + A|0|^ 
is forming, Tq signals the temperature above which there is a significant probability for 
thermal fluctuations over the central hump on the scale of the correlation length. Most 
simply, it is determined by the condition 

AViTG)ee,iTG) ^ Tg (2.2) 

where AV{T) is the difference between the central maximum and the minima of V{(f), T). 
We find \1-Tg/Tc\ = 0(A). 
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Whereas, above Tq there will be a population of 'domains', fluctuating in and out of 
existence, at temperatures below Tq fluctuations from one minimum to the other become 
increasingly unlikely. When this happens the correlation length is 

6 



UTg) = 0[ , , (2.3) 



where = is the natural unit of length, the Compton wavelength of the particles. 

It is temptingfl], |^ to identify ^eqiTc) with the scale at which stable domains begin 
to form. We shall see later that this is incorrect, for quenches that are not too slow. 
However, some care is needed if (as can happen in condensed matter physics) we never 
leave the Ginzburg regime. 

The formation of large domains is an issue that requires more than equilibrium 
physics. The most simple dynamical arguments can be understood in terms of causality. 



2.2 Causality 

We have already mentioned that causality puts an upper bound on domain size. Specif- 
ically, if G{r,t) is the two-field correlation function at time t for separation r, then G 
vanishes for r >2t approximately. This was used by Kibble to put an upper bound on 
monopole density in the early universe. If this causal bound and the Ginzburg criteria 
attempt to set scales once the critical temperature has been passed, the causal arguments 
considered now attempt to set scales before it is reached. 

Here we attempt a lower bound on domain size, an upper bound on defect density. 
Suppose the temperature T(t) varies sufficiently slowly with time t that it makes sense 
to replace V{(f), T) by V{(f), T{t)). With m?{T{t)) vanishing at T = Tc, which we suppose 
happens at t = 0, the equilibrium correlation length of the field fluctuations ^eq{T{t)) = 
\m~^{T{t)) \ diverges at T{t) = Tc. It is sufficient to adopt a mean-field approximation in 
which m'^{T) oc (T — Tc). The true correlation length ^{t) cannot diverge like C,eq{T{t)), 
since it can only grow so far in a finite time. 

Initially, for t < 0, when we are far from the transition, we again assume effective 
equilibrium, and the field correlation length ,^(t) tracks S,eqiT{t)) approximately. However, 
as we get closer to the transition ^eq{T{t)) begins to increase arbitrarily fast. As a crude 
upper bound, the true correlation length fails to keep up with C,eq{T{t)) by the time —i 
at which ^^q is growing at the speed of light, d^eq(T {—t)) / dt = 1. It was suggested by 
Kibble0 that, once we have reached this time ^(t) freezes in, remaining approximately 
constant until the time t ~ +i after the transition when is once again becomes comparable 
to the now decreasing value of ^eq- The correlation length ^eqij) = ^eq{—t) is argued to 
provide the scale for the minimum domain size after the transition. 

Specifically, if we assume a time-dependence m?'{t) = —M'^t/tq in the vicinity of 
t = 0, when the transition begins to be effected, then the causality condition gives 
tc = tQ/^(2M)-2/3. As a result, 

M^cgii) = {Mrof'\ (2.4) 
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which, with condensed matter in mind, we write as 

^=Ut)=^o[^y' (2.5) 

where Tq = = are the natural time and distance scales. In contrast to Eq.p.3|, 
Eq.|275| depends explicitly on the quench rate, as we would expect. 



2.3 QFT or Condensed Matter 

This approach of Kibble was one of the motivations for a similar analysis by Zurek[Q 
of transitions in condensed matter. Qualitatively, neither the Ginzburg thermal fluc- 
tuations, with fluctuation length Eq.pTSj nor the simple causal argument above depend 
critically on the fact that the free energy Eq. |2.1| is originally assumed to be derived from 
a relativistic quantum field theory. After rescaling, F could equally well be the Ginzburg- 
Landau free energy for the complex order-parameter field whose magnitude determines 
the superfiuid density. That is, 

F{T)= Jd'x (|^|V0r + a(r)|0p + i/?|0r), (2.6) 

in which a{T) oc irP'{T) vanishes at the critical temperature T^. The only difference 
is that, in the causal argument, the speed of light should be replaced by the speed of 
(second) sound, with different critical index. 

Explicitly, let us assume the mean-field result aiT) = aoe(T), where e = (T/Tc — 1), 
remains valid as T/Tc varies with time t. In particular, we first take a{t) = a{T{t)) = 
—aot/TQ in the vicinity of Tc. The fundamental length scale is given from Eq. |2.6| as 
^0 — ^^/2mao- The Gross-Pitaevski theory suggests a natural time-scale tq = h/ao. 
When, later, we adopt the time-dependent Landau- Ginzburg (TDLG) theory we find 
this still to be true, empirically, at order-of-magnitude level, and we keep it. 

It follows that the equilibrium correlation length Ceqit) and the relaxation time r(t) 
diverge when t vanishes as 



ee,(t) = ec 



t 



-1/2 



r(t) = To 



t 



(2.7) 



The speed of sound is c(t) = ^eq(t)/T{t), slowing down as we approach the transition as 
The causal counterpart to d^eq{t)/dt = 1 for the relativistic field is d^eq{t)/dt = 
c(t). This is satisfied at t = — t, where i = ^tqTo, with corresponding correlation length 

1/4 

^0 

(cf. Eq. ^.5| ). A variant of this argument that gives essentially the same results is obtained 
by comparing the quench rate directly to the relaxation rate of the field fluctuations. We 
stress that, yet again, the assumption is that the length scale that determines the initial 
correlation length of the field freezes in before the transition begins. 
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3 Experiments 



The end result of the simple causality arguments is that, both for QFT and condensed 
matter, when the field begins to order itself its correlation length has the form 

^ = Ut)=^o[^y- (3.9) 

for appropriate 7. 

The jump that Kibble made0 in QFT was to assume that the correlation length 
Eg . 13.91 also sets the scale for the typical minimum intervortex distance. That is, the 
initial vortex density n^e/ isQ is assumed to be 

11 1 /To \ 27 

for 7 = 1/3 and / = 0(1). We stress that this assumption is independent of the argument 
that lead to Eq.p3|. Since ^0 also measures cold vortex thickness, tq ^ tq corresponds 
to a measurably large number of widely separated vortices. 

Even if cosmic strings were produced in so simple a way in the very early universe it 



is not possible to compare the density Eq. p.lO] with experiment, in large part because of 



our uncertainty as to what is the appropriate theory. It was Zurek who first suggested 
that this causal argument for defect density be tested in condensed matter systems. 

3.1 Superfluid helium 

Vortex lines in both superfluid ^He and ^He are good analogues of global cosmic strings. 
In "^He the bose superfluid is characterised by a complex field 0, whose squared modulus 
|0p is the superfluid density. The Landau-Ginzburg theory for '^He has, as its free energy, 
F{T) of Eq. p.6| . The static classical field equation 6F/S(f) = has vortex solutions as 
before. Specifically, a simple (winding number unity) static straight vortex along the 
z-axis has the form 

$(x) = h{r)e^'', (3.11) 

where 6 = arctan(?//x) and = + y"^. For small r, h{r) = 0{r), and for large r, 
h{r) = ry(l — OI^q/t"^)), with effective width ^o- 

The situation is more complicated, but more interesting, for ^He. The reason is that 
^He is a fermion. Thus the mechanism for superfluidity is very different from that of '^He. 
Somewhat as in a BCS superconductor, these fermions form the counterpart to Cooper 
pairs. However, whereas the (electron) Cooper pairs in a superconductor form a state, 
the ^He pairs form a state. The order parameter A^i is a complex 3x3 matrix A^i- 



■^In fact, the powers of Eq j2.5| and Eq.2.8 are mean-field results, changed on implementing the renor- 
malisation group. 

^Equivalently, the length of vortices in a box volume v is 0(ndefv)- 
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There are two distinct superfluid phases, depending on how the 50(3) x 50(3) x f/(l) 
symmetry is broken. If the normal fluid is cooled at low pressures, it makes a transition 
to the ^He — B phase. 

The Landau-Ginzburg free energy is, necessarily, more complicated, permitting many 
types of vortexQ, but the effective potential V{Aai, T) has the diagonal form^] V{A, T) = 
a{T)\Aai\'^ + 0{A^) for small fluctuations, and this is all that we need for the production 
of vortices at very early times. Thus the Zurek analysis leads to the prediction Eq. p.lOl, 
as before, for appropriate 7. 



3.2 Counting vortices 



Although ^He is more complicated to work with, the experiments to check Eq. |3.10| are 
cleaner, since even individual vortices can be detected by magnetic resonance. Second, 
because vortex width is many atomic spacings the Landau-Ginzburg theory is good (7 = 
1/4). 

So far, experiments have been of two types. In the Helsinki experiment ||10|| superfluid 
^He in a rotating cryostat is bombarded by slow neutrons. Each neutron entering the 
chamber releases 760 keV, via the reaction n+^He He + 760keV . The energy goes 

into the kinetic energy of the proton and triton, and is dissipated by ionisation, heating 
a region of the sample above its transition temperature. The heated region then cools 
back through the transition temperature, creating vortices. Vortices above a critical size 
grow and migrate to the centre of the apparatus, where they are counted by an NMR 
absorption measurement. The quench is very fast, with tq/tq = O(IO^). Agreement with 
Eq. [3.101 and Eq.^]8| is good, at the level of an order of magnitude. This is even though 
it is now argued that the Helsinki experiment should no^ show agreement because of 



the geometry of the heating event. 

The second type of experiment has been performed at Grenoble and Lancaster 
Rather than count individual vortices, the experiment detects the total energy going 
into vortex formation. As before, ^ife is irradiated by neutrons. After each absorption 
the energy released in the form of quasiparticles is measured, and found to be less than 
the total 760 keV. This missing energy is assumed to have been expended on vortex 
production. Again, agreement with Zurek's prediction Eq. |3.10| and Eq. p.8| is good. 



The experiments in ^i7e, conducted at Lancaster, follow Zurek's original suggestion. 
The idea is to expand a sample of normal fluid helium so that it becomes superfluid at 
essentially constant temperature. That is, we change 1 — T/T^ from negative to positive 
by reducing the pressure and increasing Tg. As the system goes into the superfluid phase 
a tangle of vortices is formed, because of the random distribution of fleld phases. The 
vortices are detected by scattering second sound off them, and its attenuation gives a 
good measure of vortex density. A mechanical quench is slow, with tq some tens of 
milliseconds, and tq/tq = O(10^°)|. Two experiments have been performed |ll3i In 



^For '^He mean-field theory is poor, and a better value for 7 is 7 = 1/3. 
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the first fair agreement was found witli tlie prediction Eq. |3.10| , but tlie second experiment 
has failed to see any vortices whatsoever. 

There is certainly no agreement, in this or any other experiment on ^He^ with the 
thermal fluctuation density that would be based on Eq. p.3| . 

4 The Kibble- Zurek picture for the value of f is cor- 
rect. 

To do better than the simple causality arguments we need a concrete model for the 
dynamics. 

4.1 Condensed matter: the TDLG equation at early times 

We assume that the dynamics of the transition can be derived from the explicitly time- 
dependent Landau-Ginzburg free energy 

m = jd?x [^\v<p\'+am\' + \m')- (4.12) 

obtained from Eq.^]6|on identifying a{t) = a(T{t)) = aoe(t), where e = (T/Tc — 1). In a 
quench in which or T changes it is convenient to shift the origin in time, to write e{t) 
as 

e{t) = eo - —e{t) (4.13) 

for —oo < t < rQ(l + eo), after which e(t) = — 1. eg measures the original relative 
temperature and tq defines the quench rate. The quench begins at time t = and 
the transition from the normal to the superfiuid phase begins at time t = e^Tq. Times 
subsequent to that are defined by At = t — to- 

Motivated by Zurek's later numerical Q simulations, we adopt the time-dependent 
Landau-Ginzburg (TDLG) equation for F, on expressing as = (0i + 202)/"\/2, that 

-T—+Va, (4.14) 



r dt (5 „ 

where rja is Gaussian thermal noise, satisfying 

(r/,(x, t)r],{y', t')) = 25,,T(t)r5(x - y)6{t - t'). (4.15) 

This is a crude approximation for ^iYe, and a simplified form of a realistic description of 
'^He but it is not a useful description of QFT, as it stands. 

It is relatively simple to determine the validity of Zurek's argument since it assumes 
that freezing in of field fluctuations occurs just before symmetry breaking begins. At that 
time the effective potential y(0, T) is still roughly quadratic and we can see later that. 
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for the relevant time-interval —t<At<t the self-interaction term can be neglected 
(/3 = 0). 



In space, time and temperature units in which = ''"o = = 1, Eg. 4. 14 then 
becomes 

0,(x, t) = -[-V' + e(t)]0a(x, t) + 77„(x, t). (4.16) 

where f] is the renormalised noise. The solution of the 'free'-field linear equation is 
straightforward, giving a Gaussian equal-time correlation function [|1^, ^ 

{Mr,t)UO,t)) = 6abG{r,t) (4.17) 

where 

G{r,t)= dTf{t-T/2)(- — e"" e" Jo . (4.18) 
Jo KAtttJ 

and T is the renormalised temperature. At time to = coTo, when the transition begins, a 
saddle-point calculation shows that, provided the quench is not too fast. 



G{r,to)^^e-'^^^/^^''\ (4.19) 

where a = 0(1), confirming Zurek's result. 

Zurek's prediction is robust, since further calculation shows that ^{t) does not vary 
strongly in the interval —i<At< i, where At = t — to- 

4.2 QFT: Closed time-path ensemble averaging at early times 

For QFT the situation is rather different. In the previous section, instead of working with 
the TDLG equation, we could have worked with the equivalent Fokker-Planck equation 
for the probability "^[^j that, at time t > 0, the measurement of will give the function 
$(x). Thus G{r,t) of Eq. [4.17| can be written as 



6abG{r,t) = {Mr,t)UO,t)) = J P$pr[$]$„(r)$,(0). (4.20) 

When solving the dynamical equations for a hot quantum field it is convenient to work 
with probabilities from the start. Taking t = as our starting time for the evolution of 
the complex field (p suppose that, at this time, the system is in a pure state, in which the 
measurement of (p would give $o(x)- That is:- 

0(t = 0,x)|$o,t = 0) = $o|$o,t = 0). (4.21) 

The probability pt[^] that, at time t > 0, the measurement of will give the function 
$(x), is the double path integral 



pM = r^'^ * exp|z(^i[0+] - St[ct>-])\. (4.22) 

j0±(O)=*o IV /J 



10 



where 'D(j)± = Ti(j)±^iD(j)±^2) and St[4>\ is the (time-dependent) action obtained from 
Eg . 12 .11 , on substituting m{t) = m{T{t)) for m(T). 

Pt[$] can be written in the closed time-path form in which, instead of separately 
integrating (j)± along the time paths < t < tf, the integral can be interpreted as time- 
ordering of a field along the closed path C+ © C_ of Fig.l, where = 0+ on C+ and 
= 0_ on C„. When we extend the contour from tf to t = oo either 0+ or 0_ is an 
equally good candidate for the physical field, but we choose 0+. 

The choice of a pure state at time t = is too simple to be of any use. For simplicity, 
we assume that $ is Boltzmann distributed at time t = at an effective temperature of 
To = Po'^ according to the Hamiltonian corresponding to the free-field action S'[0], 
obtained by setting A = in Eq. |2.1| , in which is taken to be periodic in imaginary time 
with period /3o. We now have the explicit form for 



V(Pe 



written as the time ordering of a single field along the contour C = C-i 
extended to include a third imaginary leg, where takes the values 0+, <; 
C+, C_ and C3 respectively, for which Sc is 5'[0+], 5'[0_] and 5*0 [03]. 



(4.23) 
and 03 on 



to 



to - iPo 



Figure 1 : The closed timepath contour C_|_ © C_ , with the Boltzmann imaginary leg 

Just as we had no need to calculate explicitly in condensed matter we can 

average in QFT without having to calculate pt[^] explicitly. Specifically, 

Gab{r,t) = (<l>,(r)$,(0))t = J P$p4$]$,(r)$6(0) (4.24) 

is given by 

Gab{r,t) = (0a(r,t)0b(O,t)), (4.25) 
the equal-time thermal Wightman function with the given thermal boundary conditions. 
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Fortunately, as for the condensed matter case, the interval —t < At < i occurs in 
the /m ear regime, when the self-interactions are unimportant. The relevant equation for 
constructing the correlation functions of this one-field system is now the second-order 
equation 

for F of Eq. p^.l| . This is solvable in terms of the mode functions Xki^)^ identical for 
a = 1,2, satisfying 



Xtit) = 0, (4.27) 



subject to Xk{t) = e^*'^'"* at t < 0, for incident frequency Um = + cqM^ and 
m^(t) = e(t)M^, where e{t) is parameterised as for the TDLG equation above. This 
corresponds to a temperature quench from an initial state of thermal equilibrium at 
temperature Tq > Tc, where (Tq/Tc — 1) = cq. The diagonal correlation function G{r,t) 
of Eg . [4. 17] is given as the equal-time propagator 



G{r,t) = j fke^^-xUt)x,it)C{k) (4.28) 
^ ' dke'^xUt)x-At)C{k), 



27r^ J kr 

where C{k) = coth.{uJin{k) /2TQ)/2uJin{k) encodes the initial conditions. 

An exact solution can be given in terms of Airy functions. Dimensional analysis 
shows that, on ignoring the k-dependence of C{k), appropriate for large r (or small k), 
ieqii) of Eq.[275| again sets the scale of the equal-time correlation function. Specifically, 

G{rM)^ I dK'^^^^^F{K), (4.29) 

where -F(O) = 1 and F{k) ~ for large k. Kibble's insight is correct, at least for this 
case of a single (uncoupled) field. 



5 Vortex densities do not determine correlation lengths 
directly 

We have seen that there is no reason to disbelieve the causal arguments of Kibble for 
QFT and Zurek for condensed matter as to the field correlation length at the time of the 
transition. The excellent agreement with the ^He experiments also shows that, despite 
the very interesting simulations of Kopnin et al. ITTI , this length does, indeed, characterise 



vortex separation for condensed matter at the time when the defects form. 

However, the recent Lancaster experiment shows that this cannot always be the case. 
Significantly, for ^He the Ginzburg regime is extremely narrow, whereas for ^He it is 
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very broad. In fact, the ^He experiments begin and end in the Ginzburg regime, where 
thermal fluctuations dominate. The causality arguments are too simple to accommodate 
these facts. 

If these differences are to be visible in the formalism, it can only be through the way 
in which we relate vortex density to correlation length. We have already observed that 
the TDLG equation can be recast as the Fokker-Planck equation, whereby the ensemble 
averages can be understood as averaging with respect to the probability ]9t[$(x)] that, 
at time t, the field takes value $(x). We can use these probabilities, implicit in the 
correlation functions, to estimate defect densities. 



5.1 Classical defects in condensed matter 



It would be foolish to estimate the probability of finding profiles like $(x) of Eq. |3.11 



directly. One way is to work through line zeroes, since vortices have line zeroes of the 
complex field (f) at the centre of their cores. The converse is not true since zeroes occur 
on all scales. However, a starting-point for counting vortices in superfiuids is to count 
line zeroes of an appropriately coarse-grained field, in which structure on a scale smaller 



than the classical vortex size, is not present |I8|. That is, we do not want to entertain 
vortices within vortices. For the moment, we put in a cutoff / = 0(^o) by hand into the 
Fourier transform G{k,t) of G{r,t), as 

Gi{r,t) = J fke'^-''G{k,t)e-^'^\ (5.30) 

We stress that the long-distance correlation length C,eq(t) depends essentially on the po- 
sition of the nearest singularity of G{k,t) in the complex k-plane, independent of I. 

This is not the case for the line-zero density n^ero- For example, in our Gaussian 
approximation of the previous section Uzero can be calculated exactly from the two-point 
correlation function G{r,t) with pt[$] implicit. It can be shown, quite easily |T^, pOf that 
it depends on the short- distance behaviour of Gi{r,t) as 

the ratio of fourth to second moments of G{k,t) e~^^^^ . 

There are several prerequisites before line zeroes can be identified with vortex cores, 
and n^eroit) with ndef{t). 

• The field, on average must have achieved its symmetry-broken ground-state equi- 
librium value 

(101') = «o//?- (5.32) 

This, in itself, is sufficient to show that the causal time i is not the time to begin 
looking for defects, since (|0p) is small at this time. This, in turn, requires that 
G{k,t) be non-perturbatively (in (3) large. 
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• Only when duzero/dl is small in comparison to Uzero/^ at I = will the line-zeroes 
have the non-fractal nature of classical defects on small-scales, although vortices 
may behave like random walks on larger scales. As the power in the long wavelength 
modes increases the 'Bragg' peak develops in k'^G{k,t), moving in towards k = 0. 
This condition then becomes the condition that the peak dominates its tail. 

• The energy in field gradients should be commensurate with the energy in classical 
vortices with same density as that of line zeroes. 

We stress that these are necessary, but not sufficient, conditions for classical vortices. 
In particular, although they can be satisfied in the self-consistent linear approximation 
that will be outlined below, only the full nonlinearity of the system can establish classical 
profiles. We will term such zeroes as satisfy these conditions proto-vortices. It has to 
be said that most (but not all[^, ^) numerical lattice simulations cannot distinguish 



between proto-vortices and classical vortices. 
5.2 TDLG condensed matter 

We begin with condensed matter, which we will find to be easier. As the system evolves 



away from the transition time, the free equation Eq. [4.16| ceases to be valid, to be replaced 
by the full equation 

4(x,t) = -[-V' + e(t)+/3|0(x,t)|2]0,(x,t)+r/,(x,t), (5.33) 

where (3 is the rescaled coupling. 

In order to retain some analytic understanding of the way that the density is such 
an ideal quantity to make predictions for, we adopt the approximation of preserving 
Gaussian fiuctuations by linearising the self-interaction as 

4(x, t) = -[-V' + ee//(t)]0a(x, t) + r7,(x, t), (5.34) 

where ee// contains a (self-consistent) term O(/3(|0p)). Additive renormalisation is nec- 
essary, so that ee// ~ e, as given earlier, for t < to- 

Self-consistent linearisation is the standard approximation in non-equilibrium QFT|]35 



36| , but is not strictly necessary here, since numerical simulations that identify line zeroes 
of the field can be made that use the full self-interaction 0. However, there are none that 
address our particular problems exactly. Given the similarities with the QFT case, for 
which it is difficult to do much better than a Gaussian, there are virtues in comparing 
the Gaussian approximation for the two cases. 

The solution for G{r,t) is a straightforward generalisation of Eq. ^.l8| , 



POO _ / 1 \ '^/^ o r'T 

G(r, t)= dr fit - r/2) -— ] e"" e" Jo '^^ ^e//(t-s/2) ^ 
Jo \Atxt J 

where T is the rescaled temperature, as before. 
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Assuming a single zero of €eff{t) at t = to, a.t r = the exponential in the integrand 
peaks at r = f = 2(t — to), the counterpart of the Bragg peak in proper-time. Expanding 



about r to quadratic order gives [jT6| 



G,(0..)«f.e^/.. ■ (5.36) 

where we have put in the momentum cutoff k^^ > I = l^o = 0(.^o)of Eq. |5.30| by hand. 
For times t > e^rq we see that, as the unfreezing occurs, long wavelength modes with 
fc^ < t/rg — eo grow exponentially. 

The effect of the back-reaction is to stop the growth of ^^(0, t) — (^((O, to) = ~ 
(|0p)o at its symmetry-broken value (3"^ in our dimensionless units. A necessary condi- 
tion for this is Ymiu-^oo^effiu) = 0. That is, we must choose eeff{t) = e(t) + /3{Gi{0,t) — 
Gi{0,to)), thereby preserving Goldstone's theorem. 

At t = to; when the approximation Eq. p.3(j| is good, both numerator and denominator 
are dominated by the short wavelength fluctuations at small r. Even though the field is 
correlated over a distance ^ ^ I the density of line zeroes Uzero = 0(/~^) depends entirely 
on the scale at which we look. In no way would we wish to identify these line zeroes with 
prototype vortices. However, as time passes the peak of the exponential grows and n^ero 
becomes increasingly insensitive to /. How much time we have depends on the magnitude 
of /?, since once G{0, t) has reached this value it stops growing. The time t* at which this 
happens can be estimated by substituting e{u) for eeffiu) in the expression for Gi{Q,t) 
above. 

For t > t* the equation for Uzeroit) is not so simple since the estimate for Gi{0,t) 
above, based on a single isolated zero of eeff(t), breaks down because of the approximate 
vanishing of €eff{t) for t > t*. A more careful analysis shows that Gi{0,t) can be written 
as 

where G{T,t) has the same peak as before at r = 2{t — to), in magnitude and position, 
but G{T,t) = 1 for r < 2(t — t*). Thus, for tq ^ tq, Gi{0,t) can be approximately 
separated as Gi{0,t) = G'f^(t) + G^^(t), where 

G^^{t)=f{t) dT/[An{T + P)f^, (5.38) 



the counterpart of the Bragg 'tail', describes the scale- dependent short wavelength ther- 
mal noise, proportional to temperature, and 

roc _ 

describes the scale- independent, temperature independent, long wavelength fluctuations. 
A similar decomposition Gf/i{0,t) = Gn^^ {t) +Gn^^{t) can be performed. In particular, 
G//^^(t)/G^^(t) = 0(ri). 
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Firstly, suppose that, for t>t*, G^^{t) > GY^{t) and Gff^^{t) > GffY^{t), as would 
be the case for a temperature quench T{t) —>■ 0. Then, with little thermal noise, we have 
widely separated line zeroes, with density Uzeroit) ~ —Gff'^(t)/27iG^^(t) and duzero/dl 
is small in comparison to Uzero/l at / = ^o- Further, once the line zeroes have straightened 
on small scales at t > t*, the Gaussian field energy, largely in field gradients, is 

F^{[ d'x\{y<i)af) = -VG"{S),t), (5.40) 
Jv 2, 

where V is the spatial volume. This matches the energy 

E ^ Vn,e/(t) (271^(0, t)) = ~VG"{Q, t) (5.41) 

possessed by a network of classical global strings with density rizero-, in the same approx- 
imation of cutting off their logarithmic tails. 

From our comments above, we identify such essentially non-fractal line-zeroes with 
prototype vortices, and n^ero with Udef- Of course, we require non-Gaussianity to create 
true classical energy profiles. Nonetheless, the Halperin-Mazenko result may be well 



approximated for a while even when the fluctuations are no longer Gaussian p2 
For times t > t* 

rizeroit)^—— — -TToJ— , (5-42) 




the solution to Vinen's equation [p3 | 



= -^^-^^ef, (5.43) 

where X2 = 47r in our approximation^. What is remarkable in this approximation is 
that the density of line zeroes uses no property of the self-mass contribution to eeffit), 
self-consistent or otherwise. 

This decay law is assumed in the analysis of the Lancaster experiments. The empirical 
value of X2 used in them is not taken from quenches, but turbulent flow experiments. 



It is suggested ||14|| that X2 ~ 0.005, a good three orders of magnitude smaller than our 
prediction above. Although the TDLG theory is not very reliable for '^He, if our estimate 
is sensible it does imply that vortices produced in a temperature quench decay much faster 
than those produced in turbulence. 

Equally importantly, we shall see that, for early time at least, thermal fluctuations are 
large in the Lancaster experiments. However, for ^He, with negligable UV contributions, 
we estimate the primordial density of vortices as 

n.ero(t*)~^^ 7^./^, (5.44) 



''Calculations for X2 ior realistic values of and tq give X2 > 47r for both He and ^He 
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in accord with the original prediction of Zurek. Because of the rapid growth of G{0,t), 
{t* — to) /t = P > 1 = 0{1). We note that the factor^ of p = 8np gives a value of 
/ = 0(10), in agreement with the empirical results of |]12| and the numerical results of 



Whereas Eq. [5.44| is appropriate for ^He, the situation for the Lancaster ^ife experi- 



ments is complex, since they are pressure quenches for which the temperature T is almost 
constant at T ^ T^. Unlike temperature quenches p5|, thermal fluctuations here remain 
at full strength]^. The necessary time- independence of G^^{t) for t > t* is achieved by 
taking eeff{u) = 0{u~^). In consequence, as t increases beyond t* the relative magnitude 
of the UV and IR contributions to Gi{0,t) remains approximately constant at its value 
at t = t*. 

Nonetheless, as long as the UV fluctuations are insignificant at t = t* the density of 
line zeroes will remain largely independent of scale. This follows if Gll^^{t*) ^ G//Y^{t*), 
since Gffi{0,t) becomes scale-independent later than Gi{0,t). In []T3[ we showed that this 
is true provided 

(rQ/ro)(l-TG/r,) <C7^^ (5.45) 

where G = 0(1) and Tq is the Ginzburg temperature. With tq/tq = 0(10'^) and 
(1 — Tg/Tc) = O(10~^^) this inequality is well satisfied for a linearised TDLG theory for 
^He derivecf^ from the full TDGL theory Q, but there is no way that it can be satisfied 
for ^He, when subjected to a slow mechanical quench, as in the Lancaster experiment, for 
which tq/tq = O(10^°), since the Ginzburg regime is so large that (1 — Tq/Tc) = 0(1). 
As far as the left hand side of Eq. |5.45| is concerned, the '^He quench is nineteen orders 



of magnitude slower than its ^He counterpart. 

When the inequality is badly violated, as with ^He for slow pressure quenches, then 
the density of zeroes Udef = 0(Z~^) after t* again depends explicitly on the scale / at which 
we look and they are not candidates for vortices. Since the whole of the quench takes 
place within the Ginzburg regime this is not implausible. However, it is possible that, 
even though the thermal noise never switches off, there is no more than a postponement 
of vortex production, since our approximations must break down at some stage. The 
best outcome is to assume that the effect of the thermal fluctuations on fractal behaviour 
is diminished, only leading to a delay in the time at which vortices finally appear. Even 
if we suppose that Udef above is a starting point for calculating the density at later 
times, albeit with a different to, thereby preserving Vinen's law, we then have the earlier 
problem of the large X2 = 0{p), which would make it almost impossible to see vortices. 

For all that, a numerical simulation that goes beyond the Gaussian approximation 
specifically tailored to the Lancaster parameters is crucial if we are to understand what 
is really happening. We hope to pursue this elsewhere. 



^An errant factor of 3 appeared in the result of [|5[ 

^The temperature quench of the latter is somewhat different from that considered here, but should 
still give the same results in this case 

^°Even for ^He, T/Tc never gets very small, and henceforth we take T = TcVa Gi{Q,t) above 
^^Ignoring the position-dependent temperature of pi[ 
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6 The Appearance of Structure in QFT 



When, in Section 5.2 we set up the closed time-path formahsm for the field probabilities 
^4"^]; our aim was the limited one of establishing the role of Kibble's causal correlation 
length ^ in Eq. [4.29| . We now appreciate, from condensed matter theory, that this does 
not, of itself, imply vortices at that separation. 



6.1 Proto- vortices in QFT 

To establish a link between the correlation function G{r, t) and vortices is even more 
problematic in QFT than for condensed matter systems. Yet again, we attempt to 



count vortices by counting line zeroes|26]. In the Gaussian approximations that we shall 



continue to adopt the expression Eq. [5.31| for rizero is equally applicable to QFT. This 



counting of zeroes is the basis of numerous numerical simulations P7|, |28| , P9|] of cosmic 
string networks built from Gaussian fiuctuations. 

The prerequisites for line zeroes in condensed matter that we posed after Eq. |5.31| still 
stand for QFT (except that = M^/A), but there are further complications peculiar 

to QFT. In particular, in QFT we need to consider the whole density matrix ($'|p(t)|$) 
rather than just the diagonal elements pt[$] = ($|p(t)|$). Classicality is understood in 
terms of 'decoherence', manifest most simply by the approximate diagonalisation of the 
reduced density matrix on coarse-graining. By this we mean the separation of the whole 
into the 'system', and its 'environment' whose degrees of freedom are integrated over, to 
give a reduced density matrix. The environment can be either other fields with which our 
scalar is interacting or even the short wavelength modes of the scalar field itself PU| , PT| . 
When interactions are taken into account this leads to quantum noise and dissipation. 

In the Gaussian approximations that we shall adopt here, with ($) = 0, integrating 
out short wavelengths with k > l^^ is just equivalent to a momentum cut-off at the 
same value. This gives neither noise nor dissipation and diagonalisation does not occur. 
Nonetheless, from our viewpoint of counting line-zeroes, fiuctuations are still present 
when I = 0{M~^) that can prevent us from identifying line-zeroes with proto-vortices, if 
the quenches are too slow. 

For all these caveats, there are other symptoms of classical behaviour once Gi{0;t) 
is non-perturbatively large. Instead of a field basis, we can work in a particle basis 
and measure the particle production as the transition proceeds. The presence of a non- 
perturbatively large peak in k'^G{k; t) at k = k^ signals a non-perturbatively large oc- 
cupation number N^^^ oc 1/A of particles at the same wavenumber /cq^^. With n^ero of 



( ^.31| ) of order k^ this shows that the long wavelength modes can now begin to be treated 
classically. From a slightly different viewpoint, the Wigner functional only peaks about 
the classical phase-space trajectory once the power is non-perturbatively large p2, 33 



More crudely, the diagonal density matrix elements are only then significantly non-zero 
for non-perturbatively large field configurations oc A"^/^, like vortices. 
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6.2 Mode growth v. fluctuations 



For early times we revert to the mode decomposition of Eq. [4.27| . The term coth{uJin/2To 



appearing in it can be approximated by 2To/y/eoM. Even though this is a tempera- 
ture quench, it shows strong similarities to the pressure quench of condensed matter, 
since both the long and short wavelength contributions to G{r,t) are scaled by the same 
temperature and we cannot switch off the latter. 

The field becomes ordered, as before, because of the exponential growth of long- 
wavelength modes, which stop growing once the field has sampled the groundstates. 
What matters is the relative weight of these modes (the 'Bragg' peak) to the fluctuating 
short wavelength modes in the decomposition Eq. [4.29| at this time, since the contribution 
of these latter is very sensitive to the cutoff /. Only if their contribution to Eq. p.lO 



IS 



small when field growth stops can a network of line-zeroes be well-defined at early times, 
let alone have the predicted density. Since the peak is non-perturbatively large this 
requires small coupling, which we assume. 

Consider a quench with e(t) as in Eq. [4.13| , in which the symmetry-breaking begins 



at relative time At = t — to = 0. For a free roll, the exponentially growing modes that 
appear when At > fj^ = tqk'^/M'^ lead to the approximate WKB solution [P^] 

G{r;At) oc—-^——(^—] e ^V^Q (6.46) 



where ^^(At) = 2^ Attq/M. The provisional freeze-in time t^ when = M^/A is 

then, for Mtq < (1/A), 

MAt, ^ (MtQ)^/3(ln(l/A))2/3 ~ M^(ln(l/A))'/^ (6.47) 

where At^ = — to- This is greater than Mt, but not by a large multiple. Comparison 
with condensed matter, for which the ratio is a few (3 — 5) suggests that we don't need 
a superweak theory [p^] . 

At this qualitative level the correlation length at t^, is 

M^^'it,) ~ (MtQ)2/3(ln(l/A))^/l (6.48) 

The effect of the other modes is larger than for the instantaneous quench, giving, at 

The error term E = 0(A^/^(MtQ)^/'^(ln(l/A))^^/'^) is due to oscillatory modes, sensitive 
to the cutoff. In mimicry of Eq. p.lO| it is helpful to rewrite Eq. |6.49| as 



(ln(l/A))-i/=^[l + £;]. (6.50) 



19 



in terms of the scales tq = = M^^. The first term in Eq. |6.50| is the Kibble estimate of 
Eg . 13.101 , the second is the small multiplying factor, that yet again shows that estimate 
can be correct, but for completely different reasons. The third term shows when it can 
be correct, since E is also a measure of the sensitivity of Udef to the scale at which it is 
measured. The condition E"^ -C 1, necessary for a proto- vortex network to be defined, is 
then guaranteed if 

{tq/to)\1-Tg/T,)<C, (6.51) 

C = 0(1), on using the relation (1 — Tq/T^ = 0(A). This is the QFT counterpart to 
Eq.|05|. 

For example, suppose that this approach is relevant to the local strings of a strong 
Type-II U{1) theory for the early universe, in which the time-temperature relationship 
tT^ = TMpi is valid, where we take T = 0(10"^) in the GUT era. If G is Newton's 
constant and /i the classical string tension then, following 0], Mtg ~ 10~^A^/^(G/i)^^/^. 
The dimensionless quantity ~ 10^® — 10~^ is the small parameter of cosmic string 
theory. A value A ~ 10~^ gives Mtg ~ (Mt^,)'^,a ~ 2, once factors of vr, etc. are taken 
into account, rather than Mtg ~ 1/A, and the density of Eq. |6.50| may be relevant. 



6.3 Backreaction in QFT 



To improve upon the free-roll result more honestly, but retain the Gaussian approximation 
for the field correlation functions, the best we can do is adopt a mean-field approximation 
along the lines of I^S], as we did for the CM systems earlier. As there, it does have 
the correct behaviour of stopping domain growth as the field spreads to the potential 
minima. As before, only the large- expansion preserves Goldstone's theorem. 

G(r; t) still has the mode decomposition of Eq. [4.29|, but the modes xt ^ow satisfy 
the equation 

4, + k2 + m2(t) + A($2(0))t 



Xfc(t)=0, 



(6.52) 



where we have taken N = 2. Because A0^ theory is not asymptotically free, particularly 
in the Hartree approximation, the renormalised A coupling shows a Landau ghost. This 
means that the theory can only be taken as a low energy effective theory. 
The end result is|35|, on making a single subtraction at t = 0, is 



which we write as 



+ + m\t) + Xj fpC{p)[x;{t)x;{t) - 1] 



xtit) = o. 



Xk{t) = 0. 



(6.53) 



(6.54) 



On keeping just the unstable modes in ($^(0))t then, as it grows, its contribution to 
( |6.53| ) weakens the instabilities, so that only longer wavelengths become unstable. At t* 
the instabilities shut off, by definition, and oscillatory behaviour ensues. Since the mode 
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with wavenumber k > stops growing at time < t*, where /i^(t^) = k^, the free- roll 
density at t* must be an overestimate. 

An approximation that improves upon the WKB approximation is 

Xk{t) ^ exp / dtn{t)) (6.55) 



when 7] = M{t^ — t) > is large, and ^\{t) = fi'^{t) — k^. On expanding the exponent in 
powers of k and retaining only the quadratic terms we recover the WKB approximation 
when /i(t) is non-zero. 

The result is that the effect of the back- reaction is to give a time-delay At to t*, 
corresponding to a decrease in the value ko{t) at which the power peaks of order 

^-0(^]. (6.56, 



t* \ln{l/X) 

The backreaction has little effect for times t < t*. For t > t* oscillatory modes take over 
the correlation function and we expect oscillations in G{k;t). 

In practice the backreaction rapidly forces /i^(t) towards zero if the coupling is not 
too small For couplings that are not too weak, this requires that we graft purely 
oscillatory long wavelength behaviour onto the non-perturbatively large exponential mode 

Xtin ^ a^exp ( dtV(t')) exp ( - ^^k') (6.57) 

The end result is a new power spectrum, obtained by superimposing oscillatory behaviour 
onto the old spectrum, frozen at time t*. As a gross oversimplification, the contribution 
from the earlier exponential modes alone can only be to contribute terms something like 



^ ' ^ Mm t*) Jm<M^ 



X 



M\m{t*)\ J\k\<M' 

n(k)-w'(k) 

COS k(t - f) + ^ ' , — sin kit - t* 

k 



(6.58) 



to G, where = M{t* - ikf^/rl^'^ and W = l/4(t* - tk). The details are almost 
irrelevant, since the density of line zeroes is independent of the normalisation, and only 
weakly dependent on the power spectrum. 

The k = mode of Eq. |6.58| encodes the simple solution Xk=o(t) = a + bt when /x^ = 0. 
As observed[^] by Boyanovsky et al. this has built into it the basic causality discussed 
by Kibble0. Specifically, for r,t oo, but r/2t constant (7^ 1), 

G(r,t) ^-Q(2t/r -1). (6.59) 
r 

It follows directly that this causality, engendered by the Goldstone particles of the self- 
consistent theory, has little effect on the density of line-zeroes that we expect to mature 
into fully classical vortices, since that is determined by the behaviour at r = 0. 
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Further, for large t the power spectrum effectively has a behaviour for small k, 
unlike the white noise that would follow from Eg . 16 .461 . It has been suggested that, 
for such a spectrum, most, if not all, of the vortices are in loops, with little or no self- 
avoiding 'infinite' string (but see [^). If there was no infinite string the evolution of 
the network could be very different p7|| from that of white noise, where approximately 



75% of the string is 'infinite' p7|. Although causality due to massless Goldstone modes is 
unrealistic, the linking of causal behaviour to the long wavelength spectrum is general. 
It has to be said that this approximation should not be taken very seriously for large t on 
different grounds, since we would expect rescattering to take place at times At = 0(1/A) 
in a way that is precluded by the Gaussian approximation. 

Returning to our original concerns, if Eq. |6.51| is not satisfied, it is difficult to imagine 
how clean vortices, or proto- vortices, can appear later without some additional ingredient. 



7 Conclusions 

We examined the Kibble /Zurek predictions for the onset of phase transitions and the 
appearance of defects (in particular, vortices or global cosmic strings) as a signal of the 
symmetry breaking. Our results are in agreement with their prediction Eq. |3.9| as to the 
magnitude of the correlation length at the time the transition truly begins, equally true 
for condensed matter and QFT. 

However, this is not simply a measure of the separation of defects at the time of their 
appearance. The time i is too early for the field to have found the true groundstates 
of the theory. We believe that time, essentially the spinodal time, is the time at which 
proto-vortices can appear, which can later evolve into the standard classical vortices of 
the theory. 

Even then, they may not appear because of thermal field fluctuations. In TDLG 
condensed matter thermal noise is proportional to temperature. If temperature is fixed, 
but not otherwise, as in the pressure quenches of ^He, this noise can inhibit the production 
of vortices, although there are other factors to be taken into account (such as their decay 
rate). On the other hand, on quenching from a high temperature in QFT there are 
always thermal fluctuations, and these can also disturb the appearance of vortices. The 
condition that thermal fluctuations are ignorable at the time that the field has achieved 
the true groundstates can be written 

{TQ/roni - Tg/T,) < C, (7.60) 

where 7 = 1 for condensed matter and 7 = 2 for QFT. C = 0(1). 

This restores the role of the Ginzburg temperature To that the simple causal argu- 
ments overlooked, but does not restore thermal fluctuations as the exclusive agent for 
vortex production, as happened in early arguments. Quenches in ^He provide the major 
example for which Eq. |7.60| is not satisfied. 

What happens at late time is unclear, although for TDLG numerical simulations can 
be performed (but have yet to address this problem exactly). On the other hand, not only 
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is the case of a single self-interacting quantum scalar field in fiat space-time a caricature of 
the early universe, but it is extremely difficult to go beyond the Gaussian approximation. 
To do better requires that we do differently. There are several possible approaches. One 
step is to take the FRW metric of the early universe seriously, whereby the dissipation 
due to the expansion of the universe can change the situation dramatically ||38|| . Other 



approaches are more explicit in their attempts to trigger decoherence explicitly, as we 
mentioned earlier. Most simply, the short wavelength parts of the field can be treated as 
an environment to be integrated over, to give a coarse-grained theory of long-wavelength 
modes acting classically in the presence of noise. However, such noise is more complicated 
than in TDLG theory, being multiplicative as well as additive, and coloured [pg|, pT| , pO| . 
This is an area to be pursued elsewhere. 

I particularly thank Glykeria Karra and Eleftheria Kavoussanaki, with whom much 
of this work was done. I also would like to thank Tom Kibble and many of my colleagues 
working in this area for fruitful discussions. This is an area with a substantial, if scattered, 
literature and I have aimed to be exemplary, rather than inclusive. I apologise to any 
authors who I have not cited explicitly. This work is the result of a network supported 
by the European Science Foundation. 
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